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Graphene in spintronics Pi has so far primarily meant spin current leads of 
high performance because the intrinsic spin-orbit coupling of its tt electrons is 
very weak {2-^. If a large spin-orbit coupling could be created by a proximity 
effect, the material could also form active elements of a spintronic device such 
as the Das-Datta spin field-effect transistor Isn, however, metal interfaces often 

n 

compromise the band dispersion of massless Dirac fermions [6|]. Our measure- 
ments show that Au intercalation at the graphene-Ni interface creates a giant 
spin-orbit splitting 100 meV) in the graphene Dirac cone up to the Fermi en- 
ergy. Photoelectron spectroscopy reveals hybridization with Au5rf states as the 
source for the giant spin-orbit splitting. An ab initio model of the system shows 
a Rashba-split dispersion with the analytically predicted gapless band topology 
around the Dirac point of graphene and indicates that a sharp graphene-Au 
interface at equilibrium distance will account for only ^10 meV spin-orbit split- 
ting. The ab initio calculations suggest an enhancement due to Au atoms that 
get closer to the graphene and do not violate the sublattice symmetry. 



Graphene shows fascinating electronic properties due to its structure consisting of two 
equivalent sublattices A and B that determine its band structure with the linear dispersion 
near conical K and points in reciprocal space. A and B are identified as a pseudospin, 
and the isospin indicating K and valleys is a conserved quantum number [7] . In addition 
to the pseudo- and isospin, also the real spin is recently being considered important in 
graphene [l-^.7-13t. The main contribution to the intrinsic spin-orbit splitting originates 
from the coupling to graphene d orbitals [lot. The resulting splitting is just of the order of 
0.01 meV [Ul] which means together with the high carrier mobilities in graphene a large spin 
coherence length with measured values reaching 1.5-2 /xm ^-4]. The theoretical predictions 
are even by an order of magnitude higher so that it is at present unclear whether these 
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values represent practical limits js]. 

We have recently measured a graphene spin-orbit splitting of Ago = 13 meV Il3|]. This 
splitting occurs in the system graphene/Ni(lll) after a Au layer with one atom thickness is 
intercalated by heating between the graphene and the Ni This Au layer serves a dual 



purpose: on the one hand, it transforms the strongly bonded p graphene monolayer on 
Ni(lll) into an electronically quasifreestanding structure. This decoupling effect is apparent 
Tom the resulting linear quasirelativistic tt bands of the graphene, i. e., the Dirac cone 
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15|. 



On the other hand, the Au interlayer produces a different type of the spin-orbit splitting 
through an extrinsic effect. Indeed, the spin-orbit coupling is known to be important at 
the surface of Au [l^ and in Au nanowires |l7], and Au(lll) is frequently considered a 
model for a Rashba effect induced by breaking up-down symmetry [igIIisI . For graphene 
with its 2-atom basis (sublattices A and B) and an additionally broken up-down symmetry, 
the band dispersion including spin-orbit coupling has been calculated in Ref. [19]. In zero 
magnetic field, the dependence of the spin polarization on the two-dimensional electron wave 
vector k|| was found. The band topology was predicted to be similar to that of the unbiased 
spinless graphene bilayer but with an additional spin texture which is tangential to the 
circular constant-energy surfaces. The effect couples spin and pseudospin in such a way that 
contributions from both A an d B sublattices can lead to a substantial spin interference in 
the photoemission process 



In the present Communication, we show by spin- and angle-resolved photoelectron spec- 
troscopy that an extrinsic spin-orbit splitting Ago of the order of 100 meV can be obtained 
in graphene near the Fermi energy, i. e., an enhancement by an order of magnitude rela- 
tive to our previous finding [13]. We will show that we can very well identify the type of 
superstructure characterising the graphene- Au interface in good agreement with the bulk 
lattice constants of graphene, Au, and Ni. We determine very clearly the spin-dependent 
hybridization at the interface and its change with the superstructure. Based on ab initio 
calculations, the giant size of the spin-orbit splitting must be understood as resulting from 
more dilute Au at a closer distance to the carbon atoms than in a Au monolayer, either as 
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adatoms above or below the graphene layer. Figure 1 shows spin- and angle-resolved photoe- 
mission measurements of graphene/Au/Ni(lll). Part (a) shows the tt band with its linear 
quasirelativistic dispersion, and crosses indicate where spin-resolved photoemission spectra 
(b-d) have been measured. The splitting between spin-up spectra I^{E) (upward triangles) 
and spin-down spectra I^{E) (downward triangles) is clearly visible and amounts to 90 meV 
at k|| = 1.65 and —105 meV at ky = —1.65 A~^. Apparently, the sign of the splitting 
reverses with the sign of ky, as expected for a Rashba effect. This giant spin splitting con- 
stitutes our central experimental result. The splitting exceeds also the room-temperature 
broadening, and together with the fact that it extends to the Fermi energy, this makes it 
directly relevant for transport applications. In fact, transport properties of graphene with 
an externally produced (Rashba-type) spin-orbit splitting have been intensively investigated 



in recent years 
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23| and important predictions for such a constellation o f grap hene have 



been published: the spin Hall effect [2l| and the quantum spin Hall effect (7l]22[| and, with 



an additional exchange interaction, the quantum anomalous Hall effect 23]. 

To establish connection between our previous 13 meV data and the new 100 meV data, 
we undertook closer inspection of the spin-resolved spectra of Fig. 1. They reveal that the 
system is inhomogeneous and the high-splitting phase (^ 100 meV) and the low-splitting 
phase (^ 10 meV) are present simultaneously. Line fits and spectral decompositions (see 
Supplementary Material, Fig. S5) show that the high-splitting phase makes up between 
one third and two thirds of the spectral weight, corresponding in the simplest picture to 
half of the sample region of about 200/xm x 200/xm that is probed simultaneously. The 
decomposition is possible due to an energy shift between the two phases with the low- 
splitting phase appearing at about 200-500 meV higher binding energy. The giant size of 
the splitting and the distribution in spectral weight between the high-splitting phase and 
the low-splitting phase ranging from 1 : 2 to 2 : 1 have been reproduced in several in situ 
preparations and at different beamlines. 

We want to unveil the electronic origin of this strongly enhanced spin splitting including 
the reason for its nonuniformity and the relation to our previous results. Probing as well 
as controlling the Au layer as the likely source of the giant Rashba splitting are challen- 
ing because of its very location underneath the graphene layer. The means of control is 
basically through the amount of Au deposited initially (before heating) onto the graphene 
which typically exceeds the subsequently intercalated amount. While graphene and Ni(lll) 
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have a lattice mismatch of only 1.2%, the one between graphene and Au(lll) is much larger 
14%). Therefore, the intercalated monolayer of Au will not reach the same atom density 
as the Ni. Probing the resulting structure is possible by scanning tunneling microscopy 
(STM) because the moire effect reveals the superstructure through the beat frequency be- 
tween the graphene lattice and the Au monolayer. Moreover, the photoemission signal 
from Au5rf states under the graphene is very well visible. Figures la and 2c show that the 
graphene electronic structure resembles that of freestanding graphene. The Au layer pre- 
vents the strong carbon-Ni hybridization. Nevertheless, the graphene vr-states show also in 
the quasifreestanding phase with the Dirac cone several carbon- Au hybridization points in 
the range of 4 to 6.5 eV binding energy in Fig. 2c. Moreover, the orange arrows in Fig. 2c 
indicate further replicas of the graphene vr-band shifted to smaller and larger values of the 
wave vector ky. The shift of ky amounts to between 1/7 and 1/9 of the distance FK. The 
STM image in Fig. 2a reveals a moire pattern of similar periodicity but it can most accu- 
rately be determined from the superstructure in the low-energy electron diffraction (LEED) 
pattern as a 9 x 9 superstructure — in good agreement with the shift of the replica bands 
in photoemission. 

Dashed lines in Fig. 2c show where we probed the hybridization points by spin-resolved 
spectra (Fig. 2d). These spectra reveal that the giant spin splitting of the vr-band smoothly 
merges into an even larger spin splitting (~ 0.6 eV) of the Anbd states. This is strong 
indication that this spin-dependent Cvr-Aurf hybridization with the heavy Au is the source 
of the giant Rashba splitting in the graphene. 

Interestingly, when more Au is initially deposited, it forms after the intercalation a slightly 
different 8x8 superstructure. The Supplementary Material (Fig. SI) shows that it is char- 
acterised by an extra hybridization with other And states in the vr-band. This hybridization 
does not enhance the giant spin-orbit splitting further (Fig. S6), which indicates already 
that the differences between the 9x9 and 8x8 superstructures do not play a role for the 
giant spin-orbit splitting. 



We want to turn to ab initio theory in order to verify the giant spin splitting. Modelling 
the 9 X 9 or 8 X 8 superstructure of graphene/ Au, which would increase the unit cell by two 
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orders of magnitude over that of graphene, is unnecessary because of the similar spin-orbit 
sphtting for the two superstructures. A p{l x 1) structure and on-top position for graphene 
on 1 monolayer (ML) Au has instead been chosen (Fig. 3a). 

The distance between the graphene and Au monolayer is a parameter of enormous in- 
fluence on the spin-orbit splitting of the graphene. The fact that also the Fermi energy 
varies characteristically with this distance allows us to connect to the experiment: The ex- 
perimental Fermi level corresponds with its slight p-doping (hole doping) to 3.3-3.4 A in 
the theoretical model, and this data is shown in Fig. 3a. 3.3 A is also the equilibrium 
distance calculated before for graphene/Au(lll) [28]. According to analytical prediction, 
the band topology is peculiar around the K point with two pairs of bands — a gapped one 
and a non-gapped one This analytical model jl9l]20| does not assume any speciflc 

surface configuration. The bands have so far been confirmed for freestanding graphene in 
a supercell geometry by density functional theory calculations for an applied field E = 4.0 
V/nm [lo| but not yet for the field E realistically replaced by an interface to Au. This 
was done in Fig. 3a which shows the spin-orbit split bands at the Dirac point at K for 
the on-top geometry. In the inset, these ab initio results (symbols) are magnified near the 
Dirac point and compared to the dispersion of the analytical model (solid line), which is 
fully confirmed. This is also important for the predicted magnetic-field dependence [lol. In 
agreement with the experiment, our ab initio calculation reveals Cvr-Aurf hybridization as 
the source of the spin-orbit splitting in the graphene. The graphene vr-band {p^ orbitals) 
hybridizes indeed with the deeper Au5rf bands of d^^ and dzx type because of their matching 
symmetry. Despite the large distance of 3.3 A the hybridization is strong which has not been 
considered in the literature in connection with graphene-noble-metal interfaces. This is best 
seen after the spin-orbit coupling is turned off leaving only gaps caused by hybridization: 
we determine hybridization gaps of widths Eg^z2 ^ 0.8 eV and Eg^y^ ^ 0.5 eV (see Fig. 
Sll). Their absolute binding energies do not compare well to the experiment revealing the 
limitations of the Au monolayer as model substrate. This situation improves largely and the 
main Cvr-Aurf hybridization moves from ^ 3 eV to ^ 4 eV below Ep when the Ni substrate 
is included in the calculation (see Fig. S12). 

The spin-orbit splitting in Fig. 3a is 9 meV near Ey but increases strongly when the 
distance between graphene and Au is reduced (see Fig. SIO for the detailed behaviour 
with the distance). On the other hand, the Dirac cone is destroyed at closer graphene- Au 
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distance, giving, e. g., for 2.5 A a band gap at K of 40 meV. The reason for this is the broken 
A-B symmetry of carbon atoms in the on-top geometry (Fig. 3a) which, remarkably, does 
not manifest itself very much at larger distances. A hollow-site geometry, instead, preserves 
the A-B symmetry in the graphene, see Fig. 3b. Consequently, an intact Dirac cone is 
obtained in the hollow-site geometry also for smaller graphene- Au distances such as 2.5 A 
in Fig. 3b. We are very well able to obtain at this arbitrary interlayer distance of 2.5 A 
a giant spin-orbit splitting of ~ 70 meV. However, such close distance costs as much as 1 
eV relative to the equilibrium separation and is thus unrealistic. While the giant spin-orbit 
splitting apparently is difficult to reproduce by density functional theory in equilibrium, 
the intact Dirac cone is not. The intact Dirac cone is presently obtained with p{l x 1) 
on-top graphene/Au and has also been found for p{l x 1) graphene/Cu(lll) where the 
on-top position is determined to be energetically favourable [28]. Relative to this p{l x 1) 
on-top geometry which implies maximum A-B symmetry breaking, a moire superstructure 
such as the 9x9 one necessarily breaks the A-B symmetry to a lesser or vanishing degree 
which is favourable for obtaining an intact Dirac cone. After having investigated other 
laterally shifted positions of the p{l x 1) overlayer and a 4 x 3 moire superstructure as 
well as corrugation in the graphene all yielding spin-orbit splittings of the order of 10 meV 
at equilibrium graphene-Au distances, the conclusion is that the giant spin-orbit splitting 
will not be accounted for by the sharp graphene-Au interface alone that our structural 
characterization by LEED and STM suggests as simplest case. This means that a model 
for a realistic splitting will have to include individual Au atoms which then can obtain a 
higher coordination to the carbon and due to the resulting attraction and shorter distance 
exert a stronger spin-orbit coupling on the graphene vr-states. As seen above, critical is at 
small graphene-Au separations the preservation of A-B symmetry which leads us once again 
to the hollow site geometry. Figure 3c shows this situation for graphene/0.25 ML Au in 
a p{2 X 2) structure with the Au in the hollow site relative to the graphene. The vertical 
distance of the Au to the graphene layer is 2.3 A which (as distinct from the Au monolayer 
of Fig. 3b) is near the equilibrium determined by our structural optimization. Figure 3c 
shows that this structure enhances the spin-orbit splitting to values between 50 and 100 
meV while keeping the characteristic band topology and the Dirac point. This demonstates 
that our experimental results are also plausible theoretically. 

The 9x9 and 8x8 structures will include in the simplest case a substantial amount of 
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Au atoms arranged in the hollow sites of graphene j2J| but this does not imply a reduced 
distance, especially if one considers that 9 x 9 is also the structure which a Au monolayer 
alone forms on Ni(lll) j29|. Therefore, the 9x9 and 8x8 superstructures are probably 
not relevant for the giant splitting and both give rise to only ^ 10 meV splitting. 

The previously measured data {isl were characterized by a smaller ^ 13 meV spin-orbit 
splitting, and this low-splitting phase is principally in agreement with our present calcu- 
lations for the full Au monolayer. The published [13] band dispersion measured along FK 
reflects the presence of a sample with structural defects. It contained substantial contribu- 



tions of rotated domains visible as characteristic FM dispersions appearing along FK (see 
Fig. S2) which is not the case in the present data even when the amount of intercalated Au is 
varied systematically from zero to more than 1 ML (see Fig. SI). The intercalation process 
under the graphene, including that of Au, is a present far from understood. The accepted 



main route is via de: 
by STM directly 



'ects in the graphene 



30| , and for large molecules this can be confirmed 



3l| . The presence of many domain boundaries facilitates the intercalation 
and apparently results in the low-splitting phase. Intercalation of Au works, however, also in 
samples which are free of defects over large distances. The temperature of intercalation and 
that of the initial graphene formation are very similar in the graphene/Au/Ni(lll) system 
so that an opening and closing of the graphene appears possible during intercalation. This 
could lead to more Au locally closer to the graphene, either as subsurface Au or as adatoms. 

The p{2 X 2) plot of Fig. 3c is useful for demonstrating that an enhancement of the spin- 
orbit coupling in graphene by sparsely distributed Au atoms is possible in an equilibrium 
structure but it does not imply that the Au arrangement possesses such an ordered structure. 
We can exclude that the Au atoms substitude carbon atoms. Their size would correspond 
to two carbon atoms and the resulting distortion would manifest itself in STM. Another 
possibility are non-intercalated residual atoms above the graphene for which we have no 
direct experimental evidence. We mention that such atoms cannot be detected by our STM 
technique because they should be highly mobile on graphene and shifted along by the STM 
tip. In order to explore a realistic model which includes the Ni substrate, the band dispersion 
for a structure of graphene sandwiched between 0.25 ML Au and 1 ML Au on top of 3 ML 
Ni has been calculated as well (Fig. Sll) and gives practically the same Dirac cone with 
giant spin splitting as in Fig. 3c. 
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In summary, we report a giant Rashba splitting in graphene in contact with Au up to 
100 meV that is caused by graphene-Au hybridization. A flat Au monolayer can account 
for only ^ 10 meV spin-orbit splitting as has been reported in Ref. {ist, whereas a structure 
including laterally more separated Au adatoms residing in hollow-site positions closer to 
graphene, gives at the equilibrium graphene-Au distance rise to the enhanced splitting of 
^so ^ 50-100 meV, a realistic Fermi level position, and an intact Dirac cone. We attribute 
simultaneous presence of both 100 meV and 10 meV splittings to the coexistence of areas 
without and with extra Au either as adatoms or immersed into the graphene. 



Methods 

Angle-resolved spectra have been measured with a hemispherical analyzer and parallel 
angular detection and a 6-axes automated manipulator. Spin- and angle-resolved photoe- 
mission has been performed with a hemispherical analyzer coupled to a Rice University 
Mott-type spin polarimeter operated at 26 kV sensitive to the in-plane spin component 
perpendicular to ky. Linearly polarized synchrotron light from the UE112-PGM undulator 
beamlines at BESSY II has been used for excitation. The unmagnetized Ni(lll) surface 
was prepared as a 15-20-ML thick fllm on W(llO), and graphene was synthesized by crack- 
ing of propene at the Ni(lll) surface held at T ^ 800 K. Because the surface reactivity 
drops drastically with the graphene covering, this procedure results in exactly one graphene 
monolayer [6]. Intercalation of Au [l^ was achieved by deposition of a nominal monolayer 
coverage of Au on the graphene followed by brief annealing at 750 K. This Au coverage was 
systematically varied in a wedge- type experiment. Overall energy (of electrons and photons) 
and angular resolution of the experiments was 80 meV and 1° with both angle and spin res- 
olution and 6 meV and 0.4° with angle resolution only. The base pressure was 1-2 x 10~^^ 
mbar. 

The calculations have been conducted in the generalised gradient approximation js^] to 
density functional theory, using the full-potential linearised augmented planewave method 
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implemented in the FLEUR code 



]. We use a plane- wave cutoff of 7.37 and muffin-tin 
radii of 0.68 A for carbon, 1.22 A for Au and 1.16 A for Ni. The potential is calculated self- 
consistently with, e. g., for the calculation of Fig. 3a, 49 kn points in the two-dimensional 



Brillouin zone. Spin-orbit interaction is treated self-consistently as described in Ref. 39]. 

For simplicity, the Au monolayer was assumed to be pseudomorphic [p{l x 1)] to the 
graphene (graphene lattice constant 2.489 A) and to the Ni(lll) substrate with the same 
area density of Ni(lll). This is also the basis for the 0.25 ML Au in the p{2 x 2) structure. 
The important Au-graphene distance has been chosen in Fig. 3a based on the Fermi-level 
position and in agreement with Ref. j28| as 3.3 A. The dependence of the Fermi level and 
spin-orbit splitting on this parameter is shown in Fig. SIO. For the calculation shown in 
Fig. 3c a structural optimization has been conducted by force minimization. 

Acknowledgement: E. I. R. acknowledges funding from the NSF Materials World Net- 
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FIG. 1: Giant spin-orbit splitting Ago of the graphene tt band, (a) Angle-resolved pho- 
toemission of graphene/Au/Ni(lll) with indication by crosses where the Dirac cone of graphene 
TT-states is probed by (b-e) spin- and angle-resolved photoemission spectra (hu = 50 eV). The spin 
splitting reverses with the sign of kn as expected for a Rashba effect. 
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FIG. 2: Geometrical structure and hybridization of the intercalated Au monolayer 
under graphene. (a) Scanning tunneling microscopy and (b) low-energy electron diffraction 
show a 9 X 9 moire-type superstructure, (c) Angle-resolved photoemission reveals hybridization 
of graphene vr-bands and Au d states. Moreover, vr-bands shifted to larger and smaller ky appear 
due to the 9x9 superstructure (orange arrows). Dashed lines in (c) show at which ky points 
spin-dependent spectra are displayed in (d) {hu = 62 eV). It is seen that the spin-split tt states 
develop directly out of a large spin-orbit splitting of AuSd states 
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FIG. 3: Interface geometry and Dirac-cone splitting, (a) Ab initio calculation showing the 
Rashba-split Dirac cone for graphene and a Au monolayer in the on-top position. The approximate 
equilibrium distance of 3.3 A leads to a spin-orbit splitting of 9 meV at the Fermi level, (b) The 
Au monolayer is laterally moved to the graphene hollow sites and can now be pressed into the 
graphene to the non-equilibrium distance of 2.5 A without breaking the Dirac cone. This leads 
to a giant spin-orbit splitting of ^ 70 meV. (c) Improved model with the Au atoms still in the 
graphene hollow sites but diluted to 0.25 ML Au in a j)(2 x 2) geometry. This lowers the equilibrium 
separation between the graphene and the Au layer to about 2.3 A and causes a Fermi-level position 
similar to the experiment and a giant (50 to 100 meV) spin-orbit splitting. 
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Graphene for spintronics: giant Rashba splitting due to hybridization with Au 
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Fig. S 1: Structural characterization of the monolayer graphene/ Au phase. 

We characterize graphene on Ni(lll) intercalated with different amounts of Au by means 
of scanning tunnehng microscopy (STM) and low-energy electron diffraction (LEED): (a,e) 
Graphene on bare Ni(lll) shows a pronounced 3-fold symmetry in STM and a clear p{l x 1) 
pattern in LEED which means that graphene is perfectly in registry to the Ni substrate, (b) 
Graphene on Ni(lll) intercalated with various submonolayer amounts of Au. Underneath 
of the graphene, the Au forms islands of various dimensions and shapes. STM scans of these 



17 



islands exhibit the p{9 x 9) moire superstucture. This superstructure is attributed to lattice 
mismatch between Au and Ni. The interatomic distances are in bulk Au 2.88 A and in bulk 
Ni 2.48 A. Since the lattices of graphene and Ni(lll) match exactly, it is not surprising 
that 1 ML Au/Ni(lll) also forms a]9(9 x 9) structure at room temperature [Jac95, Ume98]. 
Underneath the graphene, this structure is rather stable since independently of the exact 
amount of intercalated Au, these p{9 x 9) islands always co-exist with areas which resemble 
cluster superlattices (marked with green arrows). Those areas can be attributed to the 
formation of an interfacial alloy between Ni and Au under the graphene. Such surface alloy 
has been observed for Au/Ni(lll) after annealing [Ves05]. However, these clusters cover only 
a minor part of the sample surface and were not found relevant to our results. The fact that 
we do not observe large areas of alloyed Au under the graphene is understandable: It was 
shown and demonstrated for CO that the presence of another species reverses the process 
again (de-alloying) [Ves05]. (c,f) Graphene on Ni intercalated with a full Au monolayer 
(nominally 1.1 ML) demonstrates in STM and LEED a perfectly periodic moire pattern. 
A quantitative analysis is easier for the LEED than for the STM since the superstructure 
can be evaluated relative to the distance between (0,0) and (1,0) spots without the need 
for calibration. Our data reveals that the moire is due to a p{9 x 9) superstructure. (d,g) 
Further increase of the Au amount leads to a p{8 x 8) structure which is less ordered in 
STM but clearly distinguishable in LEED. When the Au is deposited as a wedge, scanning 
the sample during LEED shows a clear jump between the p{9 x 9) and the p{8 x 8). (h-k) 
Corresponding development of the hybridization of graphene with Au states from a wedge- 
type sample. The p{8 x 8) phase leads to an additional hybridization around 3 eV binding 
energy. The SARPES studies of the giant Rashba splitting presented here refer to the p{9 x 9) 
phase but the Rashba splitting appears similar in the p{8 x 8) phase. 
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Fig. S 2: Original data of Fig. SI. Intensity (top) and first derivative without (middle) and 
with marks (bottom). 
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Fig. S 3: Comparison to the data of [VarOSa]. Previously we observed a Rashba splitting 
of 13 meV for grapliene/Au/Ni(lll) [VarOSa]. Fig. S3 shows published and unpublished data of 
the measurements of [VarOSa]. Apparently, the dispersions do not look similar to any of the wedge 
sample of Fig. SI. The main reason is the appearance of extra FMdispersions (no Dirac cone) 
along rK(a,c). Similarly the measurements (b,d) along rM(which should not show a Dirac cone) 
contain contributions from FK. This is due to rotational domains. The presence of domains means 
that the graphene layer was imperfect. This does not necessarily affect the SARPES measurement 
of the p-band in the vicinity of the Dirac point but spectroscopically the hybridization gaps seen 
clearly in the present data are superimposed with spurious band dispersions. Most importantly, 
the presence of domain boundaries facilitates the intercalation of the Au since imperfections of 
the graphene are considered the main pathway for intercalation of deposited metals [Ton91]. This 
leads to a different structure of the intercalated Au. The present data in Fig. Sl(j,k) do not show 
such rotational domains any more. 
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(a) Spin- and Angle-resolved measurements (b) Fitting of ARPES peak 
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Fig. S 4: Consistency of the spin-resolved and non-spin-resolved ARPES measure- 
ments. The addition of the spin-up and the spin-down SARPES spectrum leads to a broad peak 
(dashed green) in (a). On the other hand, the ARPES data [sohd blue in (a)] is sharper due to 
higher angle resolution of the ARPES setup. Therefore, in (b) we repeated the addition of two 
peaks split by 100 meV. It is seen that also in this case no splitting appears after adding the 
spectra. 
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Fig. S 5: Contributions of the phase with low spin orbit splitting to the spin-resolved 
spectra. The spectra show a rather asymmetric sphtting suggesting a less or non-spht contribution 
at higher binding energy. Such phase would correspond to the one with ^13 meV spin-orbit splitting 
in [VarOSa]. The figure shows the attempt to fit the spectra under the most simple assumption: In 
addition to the spin-up and spin-down components (here at lower binding energy) we allow for a 
non-split component (as an approximation to the small splitting of 13 meV). Because the 4 spectra 
stem from 3 measurements with slightly different ratios between the coexisting phases, their ratios 
are allowed to vary in the fit as well. The negative wave vector -1.65 A can only be reached by 
a large change in the electron emission angle. This changes also the light polarization conditions. 
Apparently, an additional component for Ni 3d states at the Fermi energy had to be introduced in 
this case. 
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(a) p (9x9) graphene/Au/Ni(1 1 1 ) (b) p (8x8) graphene/Au/Ni(11 1 ) 
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Fig. S 6: Comparison of spin-orbit splitting of p{9 x 9) and p{8 x 8) structure. Spin- 
resolved spectra are compared from a p(9 x 9) and a predominant p(8 x 8) superstructure. The 
spin splitting is very similar. 
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Fig. S 7: Au 4f core-level spectra of the intercalation process. The spectra before and 
after intercalation do not show a core-level shift which could be used to identify the sites of Au 
atoms. After graphene growth (black) the Au (^1.5 ML) is deposited at room temperature (red) 
and subsequently intercalated by annealing (blue). The Cls shift shows the successful intercalation 
whereas the Au4f does not shift visibly. The data are for Co(OOOl) instead of Ni(lll) but otherwise 
comparable. The spectra have been normalized to equal maximum intensity. 
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Fig. S 8: Cu-intercalated grapheneNi(lll) as control sample. After intercalation of a 
monolayer Cu instead of Au we do not observe a any substantial spin splitting neither in the 
graphene vr-band nor in the Cu3d states themselves. Colors distinguish spin up and spin down. 
The geometry is comparable to Fig. 1 in the manuscript. 
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Fig. S 9: Graphene grown on SiC as control sample. Spectra for two orthogonal spin 
directions in the surface plane. No spin splitting is seen, and the upper limit for the spin splitting 
from this measurement is < 9 meV parallel to kx and < 15 meV for a splitting parallel to ky. 
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Fig. S 10: Band topology and distance dependence of doping and spin-orbit splitting. 

Left: We compare the analytically calculated Rashba-type band topology (solid lines) to the cal- 
culation of Fig. 3(a) (symbols) [^^(1 x 1) on-top geometry] in the vicinity of the Dirac point and 
find full agreement. The line color gives the chirality n and the symbol color (red and blue) the 
spin polarization. The spin-orbit splitting /S^soc amounts to 9 meV. For this structure we varied 
the graphene- Au distance. Right: The spin-orbit splitting strongly depends on the distance of the 
graphene to the heavy Au atoms. This finding is in line with the results from pure metal surfaces 
such as Au(lll) and Ag(lll) where the surface and nuclear potentials were found to contribute 
multiplicatively [PetOO, Bih06, Nag09]. The sign change from n- to p-doping between d = 3.34 
and 3.41 A is in agreement with calculations for R-(v^ x ^/3) graphene/Au(lll) which give an 
equilibrium graphene- Au distance of 3.31 A [Gio08]. In addition, the behavior of the staggered 
potential (A^^) is shown. The staggered potential is induced by the Au lattice which breaks the 
equivalence of A and B sublattices of the graphene if one sublattice is in the on-top position relative 
to the Au as shown in the geometry sketch of Fig. 3(a). This effect from the staggered potential 
causes a splitting of the Dirac cone of a few meV for large graphene- Au distances but it increases 
strongly for short distances breaking the Dirac cone. 
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Fig. S 11: Hybridization as the origin of the spin-orbit splitting in the graphene. 

Calculations without spin-orbit interaction reveal anticrossings and hybridization gaps which are 
due to orbital symmetries. (The spin-orbit splitting lowers the symmetry further and introduces 
additional anticrossings.) The model is p{l x 1) graphene on a Au monolayer in the on top 
geometry. Since the vr-band is made of Pz orbitals, the observed hybridization with and dyz 
orbitals is determined by symmetry. 
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Graphene on Au on Ni (111 ) 




Fig. S 12: Effect of the Ni substrate on the electronic structure at the graphene/Au 
interface. The spin-orbit splitting in the graphene clearly stems from the contact with the heavy 
Au atoms and we restricted the modeling in the manuscript to this bilayer of graphene and Au. 
Nevertheless, we have also investigated which principle influence the Ni substrate underneath has. 
Colors mark states with a high probability at the graphene (red), Au monolayer (violet), and Ni 
substrate (grey). The experimental lattice constant of Ni(lll) agrees within 1% to that of graphene 
and was chosen to be equal here. The Au monolayer was compressed to the same lattice, as before. 
The distance graphene-Au was set to 3.3 A. This distance nicely reproduces the position of the 
graphene states from the experiment. The same distance was used in the graphene-Au bilayer 
model of Fig. SIO but we observe that between these two models the binding energies of the Au5d 
bands suffer considerable changes due to the Ni substrate, mainly Au-Ni hybridizations which make 
the Au5d bands difficult to distinguish. The dzx band from Fig. Sll originating at -2 eV at F and 
arriving at -6 eV at K can be distinguished also here. A new Au5d state appears at -1.5 eV at K. 
The Au sp band of Fig. Sll is not observed in the experiment, and the present figure shows that 
the reason is that this band is very strongly hybridized with the Ni and disappears as a well-defined 
band. 
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Fig. S 13: Same model as in Fig. S12 but with additional Au adatoms. This calculation 
explores again the influence of the Ni substrate (grey symbols) under the Au monolayer (violet 
symbols) on the graphene vr-states (red and blue symbols). In this case, a giant Rashba splitting 
similar to Fig. 3(c) is achieved by adding p(2 x 2) 0.25 ML Au as adatoms on the surface in 
graphene hollow sites in the geometry of Fig 3c. In this way the graphene is sandwiched between 
0.25 ML and 1 ML Au. The upward and downward direction of the triangles marks the spin, in 
particular the upward red and downward blue triangles for graphene states. Left: The states due 
to the additional Au adatoms (black triangles, Au') can be seen to interact with the graphene 
states through different hybridizations than the Au monolayer. Right: The spin-orbit split Dirac 
cone from the graphene- Au bilayer model of Fig. Ic is thus conflrmed in the present, more realistic, 
conflguration. The fact that higher and lower energies are not symmetric about the Dirac point 
is not due to the Ni substrate but to the p{2 x 2) 0.25 ML Au adatoms, like in Fig. 3(c) of the 
manuscript. 
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Fig. S 14: Origin of the spin-orbit splitting in the Au. Comparison of +A:|| = 0.6 A and 

—fell data for graphene/Au/Ni(lll). The behavior of the spin polarization [p = (/^ — /^)/(/^ + /^)] 
shows an almost vanishing spin polarization at fcy = 0, as expected for a Rashba effect, and a 
rather clear reversal for — /cy if we take into account that the spin polarization of Au5d states is 
also subject to spectroscopic effects such as the linear dichroism of Au5d emission [Sto91] which 
do not reverse with the sign of /cy. This means that the Rashba-type spin-orbit splitting of the 
graphene tt states is accompanied by a Rashba effect on the Au monolayer itself. Such effect has 
similarly been observed for Au/W(110) [Shi08, VarOSb]. 
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